ON THE COMPUTATION OF SPECTRA IN FREE PROBABILITY 



FRANZ LEHNER 

Abstract. We use free probability techniques to compute borders of spectra of non her- 
mitian operators in finite von Neumann algebras which arise as "free sums" of "simple" 
operators. To this end, the resolvent is analyzed with the aid of the Haagerup inequality. 
Concrete examples coming from reduced C* -algebras of free product groups and leading 
to systems of polynomial equations illustrate the approach. 



1. Introduction 

Computation of spectra of convolution operators on discrete groups was one of the 
motivations behind the development of free probability. It has some interest for harmonic 
analysis and analysis of random walks on free products of discrete groups, where under the 
name of the transition operator these operators carry much information about the random 
walks under consideration. Free probability is an abstract framework for harmonic analysis 
on the free group in the language of non-commutative probability. Let us introduce briefly 
the terminology of the latter, which is used throughout the paper. 

Definition 1.1. A non- commutative probability space is a pair (A, (p) of a (complex) al- 
gebra A with unit / and a linear functional tp on A which satisfies ip(I) = 1. We will 
usually work with C*-algebras and faithful states, i.e. positive unital linear functionals. 
The elements of the algebra are called (non- commutative) random variables. The dis- 
tribution of a random variable a 6 A is given by the collection of its moments (p(a n ), 
n = 0, 1, 2, ... . If a is self-adjoint, then this corresponds to a probability measure on 
the spectrum of a. The distribution of a family of random variables is the collection of 
their mixed moments, which abstractly can be interpreted as a linear functional on the 
polynomials on noncommutative variables. 

Given a noncommutative probability space (A,ip), the subalgebras Ai C A are called 
freely independent (or free for short) if 

(1.1) ^( a i a 2 • • -a n ) = 

whenever aj G Ai i with ip(aj) = and ij ^ for j = 1, . . . , n — 1. 

In the rest of this section we review some of the necessary facts from free probability 
theory and refer the reader to [|VDN92|| and |PPOO|| for further information. 



One of the basic problems is free convolution: Given free non-commutative random 
variables a, b e A whose moments are known, compute the moments of their sum a + b. 
From fll.lj ) one sees immediately by induction that the mixed moments of a and b only 
depend on the individual distributions of a and b. Therefore the moments of a + b only 
depend on the individual moments of a and b, and thus the term free convolution of 
distributions is justified. In the case of selfadjoint random variables, once the moments 
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are known, one can proceed to compute the spectra, see section 0. What to do in the 
non-selfadjoint case will be the subject of the rest of this paper. 

The computational machinery of free convolution was found independently and about 
the same time by W. Woess ||Woe86|| , D. Cartwright and P. Soardi ||CS86|| , J.C. McLaughlin 



McL86|| in the language of random walks, and in most generality by D. Voiculescu ||Voi86 
Using the conventions of the latter, the recipe goes as follows. The Cauchy transform of 
a random variable a is the function 



^ oo 

;i-2) G a (o = ^((c-a)- 1 ) = -E 



S 71=0 S 

which is defined and analytic at least for \(\ > \\a\\. It has an inverse (under composition) 
in some neighbourhood of infinity which has the form 

(1.3) K a {z) = G- a l {z)= l -(l + R a (z)) 

where R a (z) = c±z + c 2 z 2 + . . . is analytic and one has 

(1.4) R a+b {z) = R a {z) + R b {z). 

This allows in principle to compute the moments of the sum a + b. However, function 
inversion is difficult and for non-selfadjoint operators knowledge of moments is not suffi- 
cient to determine the spectrum. For the latter, a more detailed analysis of the resolvent 
is necessary. 

The paper is organized as follows. 

In section |2| we review the case of self-adjoint operators, where the computation is 
reduced to the solution of a moment problem. 

In section [3] we use Haagerup inequality to estimate the norm of the resolvent of a sum 
of free operators. 

In section f| we look at the sum of two free operators and give a somewhat easier proof 
for this case. 

Finally in section |5] we present some examples of computations. 

2. Self-adjoint operators 

We collect here some well known techniques about computation of spectra of sums 
of self-adjoint free operators. The main tool to study spectra of selfadjoint operators is 
the spectral measure and its Cauchy- Stieltjes transform. Recall that to every selfadjoint 
operator T in a C*-probability space with faithful state ip one can associate a probability 
measure /i which is characterized by the property that <p{T n ) = J t n dfi(t). By faithfulness 
of the state <p one has supp fi = cr(T). Thus knowing this measure one knows a fortiori 
the spectrum of T. In order to compute this measure one has to solve a moment problem. 
This is usually done via the Cauchy transform ( pL2|) , which can be written in this case 

_ f dfj,(t) 

This is an analytic function on C \ cr(T) mapping the upper half plane to the lower half 
plane. The absolute continuous part of the measure can be found by the Stieltjes inversion 
formula (cf. |EV9§) 

dt 7T e\0 V ' 
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Atoms can be detected by studying the poles of G((), namely 

MW)= Km (C-t)G(C)- 

s 

non-tangentially 

When studying free sums of self-adjoint operators, by ( |1.4j ) one usually is given the 
inverse K(z) = -(1 + Yl^-ii 2 )) °f the Cauchy transform G(() and the latter is not 
accessible directly, but it is still possible to compute the spectrum. We are grateful to 
U. Haagerup for showing the following lemma to us. 

Lemma 2.1. Let A be a C* -probability space with faithful state tp, and let L 2 (ip) be the 
associated L 2 -space (the closure of A under the Hilbert norm \\X\\ 2 = V (X*X) 1 ' 2 ). LetT 
be a normal operator, then cr(T) = {A : ||(A — T) _1 || 2 = oo}. 

Proof. Given A G cr(T) and e > we have to find A G cr(T) such that |A — A | < e 
and || (A — T) _1 || 2 = oo. To this end consider the square Q of diameter 2e centered at 
Ao- The spectral measure of this square is positive, say ^t(Qo) = S > 0. If follows that 
one of the four quarters of Qo has mass at least 5/4. Let Q\ be such a quarter. Doing 
the same argument again, we find a subsquare of Q\ of mass at least 5/16. Repeating 
this construction gives rise to a sequence Q n C Q with diam(Qn) = e/2 n ~ 1 and mass 
UriQu) > S/4 n . Let A be the limit point of the sequence Q n , then every point r\ e Q n is 
at distance \tj — A| < e/2 n from A and thus 




oo 



□ 



It is therefore enough to know the L 2 -norm of the inverse, which in the case of a 
self-adjoint operator is easy to compute. For Im£ ^ it is 



l(C-T)- 1 " 2 



1 1 \c-t\ 2 



1/1 1 



c-cvc-t c-t 

G(C) - G(C) 



dfiTif) 



c-c 

while for Im ( = we take the limit to get 

IKc-t)- 1 !!^-^). 

In terms of the inverse function K(z) (11.3) this becomes 



„ \— -r^\ =7 s Imz^O 

K(z) — T)~ \\ - J K(z)-K(z) 



K(z)-K(z) 
1 

'K'(z) 



Im z = 



For concrete examples of the use of this see e.g. ||AK88|, |FTS94 . 
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3. The non-selfadjoint case 



Spectral theory of non-selfadjoint operators is complicated by the absence of the spectral 
measure^ and one has to analyze the resolvent directly. We will see that the L 2 -norm of 
the resolvent of free sums is easily computable and Haagerup inequality will provide a 



replacement of lemma |27[ 



3.1. The free resolvent. The resolvent for the sum of convolution operators on free 
products of discrete groups has been computed by several authors, among the first are 
|Woe86| | , | fStc85| |, [ |Soa86| and | |CS86| | . On the free group itself it was also computed in 
Aom84 ] and [GW8b'|. The resolvent of a sum of two free operators has been computed 



in ||Klaa92|| and very elegantly in [|Haa97|| , see sec. |] below. 

Another proof for the free resolvent formula and the i?-transform can be based on the 
following lemma (cf. ||FTS94|| , ||Leh99|| ). 

Lemma 3.1. Let Si,...,Sn E B(H) be arbitrary operators and assume that the sum of 
alternating products 



E E 

n=l t 1 ^i 2 =/=---=/=i n 

(the sum over all products where neighbouring factors are different) converges, then it 
equals 

N 
i=l 

Convergence criteria will be presented below. The following is a simple reformulation 
of this summation formula. 

Proposition 3.2. Let X\, . . . be operators on Hilbert space. Assume that I — Xi is 
invertible for every !£{!,..., N} with inverse (I — Xi)^ 1 = I + Si. Then 



E* 



-i 



q. q. . . . q. 



'+E E 

n=l i 1 =/=i2^--=/=t n 

provided the sum on the right hand side converges. 

Sums over alternating indices are well suited for free probability, as the very definition 
of freeness ( |1 . 1|) suggests and we have the following theorem. 



Theorem 3.3. Let (A, (p) be a non- commutative probability space an let Ti E A be freely 
independent random variables with Cauchy transforms Gi(() = (f((( — Tj) -1 ). The Cauchy 



transform is invertible at infinity with an inverse of the form K^z) - 
<f((Ki(z) — = z and thus 

{K^-TA-^zil + S^z)) 

with <f(Si(z)) = 0. Then the resolvent of T = Y^Tj at K{z) = -(1 
written formally as an infinite sum 

(oo 
n=l i^i^^—^i^ 



Ri(z)), i.e. 



Ri( z )) can be 



5, 



some cases Brown's spectral measure [Bro86] is computable and can serve a replacement - see 
| HL98| , |BL99| 
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In particular, the expectation of the resolvent is ((K(z) — = z an d ( EH) holds. 

If in addition 

(3-2) y4i<! 

then the sum ( |3.1| ) converges in L 2 (ip) and its norm is 



(Kiz)-^)- 1 9 = M 2 ( 1 + E E ll^ilBll^llS-ll-si 

n=l i 1 =£i 2 ^---=£i 
2 



2 

!n 112 



(3.3) =N a ll-y l|gill2 2 

I ^1 + 11^11^ 



-i 



-i 



IIWO-t,: 



-i 



Moreover, boundedness of\\(K(z) — l \\2 i m P^ es boundedness of\\(K(z) — J^Tj 

(operator norm). 



-H 



The reader easily verifies the following reformulation in terms of moment generating 
functions. 

Corollary 3.4. Let a%, % = 1, . . . , n be freely independent random variables in the non- 
commutative probability space (A,(f). Let fi(s) = -Gi(-) = <p((l — saj) -1 ) the moment 
generating function and denote 

Oi(s) = (1 - SG^r 1 - fi(s) 

the centered part of the resolvent. Given A G C, assume that there exist numbers z / 
and Si such that Sifi(si) = z for each i — 1, . . . , n and such that A = - + — -) 

i. if 

< 3 - 4 ' E < 1 



(7. 



* 112 



tnen A G" c(X] a « 
2. // 



II ° ||2 

E \\ a i \\2 _ i 
1 , II? II 2 ~~ 



± \ ||U/j||2 

and £nere are A' arbitrary close to A /or which ( |3.4| ) holds, then A G er(X^ a i)- 

Theorem |3l^ provides an easy to check sufficient criterion for the boundedness in L 2 ((p). 
By virtue of the following variant of the Haagerup inequality we will see that it is also 
sufficient for boundedness in operator norm. 

Proposition 3.5 ( |pHR97| , Lemma 3.4]). Let (A,r) be a tracial C* -probability space and 

o 

let Ai be free subalgebras with orthonormal bases Xj = {J, x^i, x^2, ■■■}■ Let X{ — Xj\ {/} 
be the centered part and 

o 

X = {Xi^k 1 X2,k 2 ■ ■ - x n,k„ '■ Xjfe X ip ij ^ ij+l} U {/} 

the free product of these bases. This is an orthonormal basis of the free product *Ai and 
can be decomposed X = M E n into the subsets E n of words of length n. For a finitely 
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supported operator a G spanX denote by Fi(a) its i-support, i.e., the set of all x G X{ 
appearing in the words of the expansion of a. Then for any a G spani£ n we have 



1/2 



Ml < (2n + 1) max I 

, x£Fi(a) 



Proof of Theorem \3. 4 Formula ( |3.1|) follows from Proposition |37|. For the L 2 -norm, 
note that the freeness condition implies that for different sets of indices ii,i 2 , . . . ,i m 
and ji, j 2 , . . . ,j n the summands S^S^ . . . S im and Sj-^S^ . . . Sj n are orthogonal and that 



«m 112 



I'-'ii II 2 ll^i: 



'12 112 



1 Si m || 2 - Then formula ( [3.3|) follows from Lemma |0 



H'S'ilS^ . . . S; 

To show boundedness of the resolvent in operator norm, we resort to analytic func- 
tions, as in | FTS94 , Chapter 2, Lemma 1.6]. Assume that the L 2 -norm is bounded, 
i.e., ll^nll^ • • • ll^nll^ < o°. Then we can define an analytic function on the open disk 
D 



n=l i 1 =£i 2 =£---=£i n 

which by the summation formula from Lemma |3.1| equals 



2 • • • ll^in II 2 



Fit) 



* 112 



1112 



This is a rational function and by assumption it has no pole on the circle {£ : |£| = 1}. 
As a rational function it has finitely many singularities and is therefore analytic on some 
disk of radius 1 + e with e > 0. It follows that there exists some constant C independent 
of n such that the Taylor coefficients satisfy 



?n 1 1 2 



< C(l + s)- n 



Now with the help of the Haagerup inequality (Proposition |3.5|) we obtain the estimate 



I + Sii Si 2 ■ ■ ■ Si n 

n=l i 1 ^i 2 ^—^i n 



< l + ^(2n+l)max||^|| ( ^ \\SiX> 



n=l 
00 



< 1 + ^(2n + 1) max H^H + e)^ n/2 



n=l 



< OO 



1/2 



15, 



In II 2 



□ 



The criterion for boundedness of the L 2 -norm 



(3.5) 



\Si 



i\\2 



\Si 



< 1 



HI2 
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can be reformulated in terms of symmetric functions. Recall that the elementary sym- 
metric functions in n variables X\,X2, ■ ■ ■ ,x n are defined as 

E)~ Efc {x\ j X2j ■ ■ ■ ; 3^n) ^ %ii ■ ■ ■ -^ij. j 

»i<—<t* 

their generating function is 

n n 
i=l k=l 

The left hand side of (|3.5|) is a rational symmetric function of Xi = \\Si\\ 2 and can be 
expressed in terms of the elementary symmetric functions as follows. 

U(l +tXi) 

U0- + tXi) 



2 

at t = 1 and x^ = \\Si\\ 2 we get therefore the condition 

Ylk=ik Ek{\\Si\\ 2 2 , . . . , \\s n \\l) 



1 + J2k=i E k{\\Si\\ 2 ,...,\\s n 



|2> 



< 1 



which is equivalent to 



£)(*-l)£l fc (||5 1 ||;,...,||5 B ||5)<l 



k=l 



For n = 2 this reduces to the simple condition j | aS'i ] | ^ \\S2W2 < 1- ^ e wm investigate this 
case in section f| below. 

4. Case of two operators 

In the particular case of the sum of two operators the analysis can be somewhat sim- 
plified and there are shorter proofs. We can use the following proposition instead of the 
Haagerup inequality from Proposition \5.5\ It is much easier to prove. 



Proposition 4.1 ( |[HL98| , Prop. 4.1]). Let («M,t) be a non- commutative von Neumann 
probability space with faithful trace state t and let a, b G M. be arbitrary centered (r{a) = 
r(b) = 0) *-free random variables. Then the spectral radius of their product is 

p{ab) = \\ab\\ 2 = \\a\\ 2 \\b\\ 2 



Corollary 4.2. For a and b as in Proposition \4 . i| , the circle of radius p(ab) = \\ab\\ 2 is 
part of the spectrum o~(ab). 

Proof. We can rescale a and b and assume that ||a|| 2 = ||6|| 2 = 1. We have to show that 
1 — tab is not invertible whenever \t\ = 1. whenever \t\ < 1, then 1 — tab is invertible and 
t is in the complement of the spectrum. By orthogonality, the L 2 -norm of the inverse is 
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and this grows unboundedly, as \t\ tends to one. Consequently the operator norm of the 
inverse becomes unbounded as \t\ —>■ 1. Since the resolvent is continuous (even analytic) 
on the complement of the spectrum, any number t of modulus \t\ — 1 is in the spectrum 
of ab. □ 



Part (|) of the following proposition is taken from | Haa97| . 
Proposition 4.3. Let a, b E A and \s\ < -J^, \t\ < f(s) = (f((l — sa)^ 1 ), g(t) 
<£>((! — tb)^ 1 ). Put a(s) = (1 — sa)' 1 — f{s), b(t) = (1 — tb)' 1 — g(t) and assume that 



(4.1) 
Set 

(4.2) 

then 
ft) 



sf(s) 



tg{t) + 0. 



sf(s) 



X — a — b 



m (1 



— sa) 



1 - 



a(s)b(t) 



[1-tb). 



(ii) If \\a{s)\\l \\b(t)\\l < \f{s)\ 2 \g{t)\ 2 , then X G C \ a (a + b). 
(Hi) If \\a(s)\\l \\b(t)\\l = \f(s)\ 2 \g(t)\ 2 , then X G a(a + b). 



Proof. Part (|) can be derived from Theorem |3.3| and can also be verified directly by 
expanding the right hand side. Parts (U) and ([Hi]) follow from Proposition ]4.1| and Corol- 
lary |4.2j respectively and the fact that 1 — sa and 1 — tb are invertible. □ 



5. Examples 

The most interesting examples are perhaps those coming from convolution operators 
on free products of discrete groups. Let us consider here non-self adjoint convolution 
operators supported on the generators of free products of cyclic groups. For a survey on 
spectra of such operators on more general finitely generated discrete groups see [|dlHRV93a] , 
|:llHRV93b|1 . Working with finite groups has the advantage that the involved equations 
are algebraic and the powerful machinery of algebraic geometry is available to obtain 
quite explicit results. Let us therefore recall some facts from algebraic geometry. For 
more background on algebraic equations we refer to the excellent survey ||CL098|| and 
text book ||CLU97|| . 

We used Singular [ pPS98|| for the algebraic computations, GNU octave ||Eat98|| for the 
eigenvalue computations and Mathematica ||Wol96|| for visualizations. 



5.1. Eigenvalue approach to polynomial equations. The most efficient numerical 
approach seems to be the matrix eigenvalue method following | ULU98 , section 2.4]. Recall 
that systems of polynomial equations are in one-to-one correspondence with polynomial 
ideals. Let fx, ■ ■ ■ , f m £ C[xi, . . . ,x n ] be polynomials and denote by V = V(fi, . . . , f m ) 
the set of solutions (also called algebraic set or algebraic variety) of the system of poly- 
nomial equations 



(5.1) 



/i = 0,...,/ n 



0. 
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On the other hand, on can associate to this system the ideal 

I = (fl,---,fm) = {J^ft/i : 9i £ C[x U . . . ,X n ]j 

it generates. Then it is easy to see that the variety V does not depend on the generators, 
but only on the ideal I. That is, V = V(I) = {x G C : f(x) = V/ G 1} (note however, 
that different ideals can lead to the same variety V). 

There are close relations between the Properties of the ideal and its variety. In particu- 
lar, one can show that the system ( |5.1|) has finitely many solutions (V is zero- dimensional) 
if and only if the quotient algebra A = C[xi, X2, ■ ■ ■ , x n ]/I is finite dimensional. Then for 
any polynomial h G C[xi, . . . , x n ], the eigenvalues of multiplication map 

M h : A -> A 
[g] >-> [hg] 

coincide with the values of the polynomial h(£), evaluated at £ G V. In particular, 
choosing h = x\ one can compute the coordinates of the solutions; in fact the eigenvalues 
and eigenvectors of a single matrix Mh are sufficient to determine all the solutions, see 
[ |GL098| , section 2.4]. Such an approach avoids the propagation of rounding errors which 
are common to elimination methods. 

In order to do concrete computations in the quotiont algebra A one needs a normal 
forms for its elements. This can be accomplished using Grobner bases. Given any (linear) 
order on the monomials, one can use a generalized Euclid's algorithm to write any / G 
C[xi, . . . , x n ] as a sum 

/ = H \-q m fm + r 

The remainder r of this division is a natural candidate for a normal form of / modulo the 
ideal J, however it may not be unique and there may exist / G / with r / 0. Fortunately 
there always exist generating sets which are nice in this respect. A Grobner basis of the 
ideal J is a generating set G = {gi, . . . ,g p } s.t. the remainder r of Euclidean division is 
unique, i.e., / G I if and only if r = 0. Grobner bases always exist and can be constructed 
with Buchberger's algorithm. With such bases the remainder r can serve as a normal form 
for the elements of A. Another characterizing property of Grobner bases is the following. 
For a polynomial / G C[x%, . . . , x n ] denote L(f) the leading monomial of / (with respect to 
the chosen monomial order) and let L(I) be the ideal generated by the leading monomials 
of all elements of I. Then G is a Grobner basis if and only if (L(G)) = L(I). In the case 
of a zero dimensional ideal, a natural basis of A = C[x\, . . . ,x n ]/I are the equivalence 
classes 

B = {[x a ] : x a (L(G))} 

i.e. the (finitely many) equivalence classes of monomials which are reduced with respect 
to the Grobner basis G. 



5.2. Elimination. Grobner bases are also needed for elimination. Roughly speaking, 
elimination of variables corresponds to projection of algebraic varieties. Elimination is 
done by computing the intersection I k — I fl C[xk+i, ■ ■ ■ , x n }. For suitable monomial 
orders (e.g. lexicographical order), one can use the fact that a Grobner basis G of I 
has the property that Gk = G fl C[xk+i, ■ ■ ■ ,x n ] is a Grobner basis for J^. However, 
lexicographical Grobner bases are expensive to compute and it is usually more efficient to 
work with other monomial orders, like degree reverse lexicographical order and use other 
methods for elimination. 
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5.3. Free product groups. For simplicity, let us consider the sum of two convolution 
operators u m + v n , where u m , v n are the generators of C|(Z m * Z n ). The moments of u m 
are 

1 k = mod m 
otherwise 



and the inverse of 1 — su m is 

x _ 1 + su m + • • • + s^uZ' 1 



1 — su 



1 - s m 



Thus for the moment generating function we obtain 

1 



f(s) = r((l-su m )- 1 ) 



1 - s rn 

From Proposition [4. 3| we infer that if the equations (|4. 1|) and (|4.2|) satisfied, which in our 
particular example read 

s t 



1 - s m 1 - t n 
A = - + s™" 1 

and if s and t are such that the L 2 -norm condition 



\s 


2 + 


s 


*+... 


+ 


s 


2m-2 


t 


2 + 


t 


* + ••• 


+ 




2n-2 




1 - s m 


2 


i - t n 


2 



< 



1 - s m \ 2 II -t n \ 2 



holds, then A is not in the spectrum of u m + u n . In order to find the outer border of the 
spectrum (that is, the curve where the L 2 -norm becomes infinite), we have to solve the 
equations 

(5.2a) s(l - t n ) - t(l - s m ) = 

(5.2b) At - 1 - s m ~H = 

(5.2c) (\s\ 2 + \s\ 4 + ■ ■ ■ + \s\ 2m - 2 ) (\t\ 2 + \t\ 4 + ■ ■ ■ + \t\ 2n - 2 ) < 1 

and look where the inequality (|5.2cj ) becomes an equality. There are two approaches to 



cope with the non- algebraic part ( 5.2c|) . One is to solve the algebraic part ( 5.2a ), ( 5.2b ) 



first and verify (|5 . 2 cj ) numerically. The alternative approach is to introduce new variables 
for the real and imaginary parts of A = x + iy, s = si + is2 and t = t\ + it 2 , and to 
separate real and imaginary parts of the equations, which now become purely algebraic, 
but with more unknowns. The real solutions of the new system correspond to the complex 
solutions of the original one. 

The second approach makes sense if the larger system is sufficiently small to allow 
complete algebraic elimination of the parameters Sj and tf, for a numerical solution the 
first approach is preferable. 

Example 5.1 ("«2 + v^'). Let us consider the simplest non-trivial example (m = 2 and 
n = 3) and compute the border of the spectrum of 112 + v%. If for a given A the system 

(5.3a) s(l -t 3 ) -t(l - s 2 ) = 

(5.3b) At - 1 - st = 

(5.3c) |s| 2 (|t| 2 + |t| 4 ) < 1 

has a solution, then A is not in the spectrum of u 2 + v$. 
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Figure 1. Spectrum of u 2 + «3 



In order to get the border of the spectrum, we replace the inequality ( |5.3c|) by an 
equality and separate real and imaginary part, which leads to the system 



(5-4) 

— S\ + 3 S2 ti t 2 + 3 S\ t\ t 2 — ^2 ^2 ~l~ S l ^1 — S 2 ^1 — ^ Si $2 ^2 "I" s l — ^1 = 

— s 2 t\ — 3 Si t 2 t 2 "I" 3 S2 t\ ^ + Si t 3 + 2 Si S2 tl + S 2 ^2 — S 2 ^2 "I" s 2 — ^2 = 

— x ti + si ti + y t 2 — s 2 1 2 + 1 = 
—y h + s 2 U - x t 2 + si t 2 = 

S-^ — \~ S<2 t-^ ~\~ 2 S t-^ is<2 ~\~ 2 S<2 ij<2 ~\~ S-^ ij<2 ~\~ ^2 ^2 ~~ ^~ ^1 ^1 ~~ ^2 ^1 ~~ ^1 ^2 ~~ ^2 ^2 — 

Using degrevlex order there is a Grobner basis of size 34 for the ideal I generated by 
the left hand sides of these equations and elimination of Sj and t» succeeds: the ideal 
JnC[i, y] is generated by a single polynomial and we get the implicit equation 

x 16 + 8 x 14 y 2 + 28 x 12 y 4 + 56 x 10 y 6 + 70 x 8 y 8 + 56 x 6 y 10 + 28 x 4 y 12 + 8 x 2 y 14 + y 16 

-6 x 15 + 54 x 13 y 2 + 226 x 11 y 4 + 238 x 9 y 6 - 18 x 7 y 8 - 158 x 5 y 10 - 74 x 3 y 12 - 6 x y 14 



.14 



+5 x 
+32 x 



65 x 12 y 2 + 305 x 10 y 4 + 435 x 8 - 225 x 6 y 8 - 179 x 4 y 10 + 107 x 2 y 12 + y 14 
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400 x 11 y 2 - 400 x M y 4 + 480 x' y° + 384 x° y 8 - 80 x A y w -Wxy 

12 i r 7Am lQ rii 2 1nQQ J„,4 K/1Q^6„,6 , rA*7„4„,8 7 n „2 „,10 , „,12 



,5 „,8 



„3 „,10 



J. 2 



-59 x iZ + 74 x iU y z - 1033 x 8 y 4 - 548 x 6 + 547 x 4 y* - 70 x z y w + y 



-40 x 11 + 776 x 9 y 2 - 1008 x 7 y 4 + 112 x 5 y 6 - 104 x 3 y 8 + 8 x 



10 



+136 x 10 + 48 x 8 y z + 736 x° y" - 176 x 4 + 24 x z y 6 

-32 x 9 + 224 x 7 y 2 - 224 x 5 y 4 + 32 x 3 y 6 

-48 x 8 - 32 x 6 y 2 + 16 x 4 y 4 = 0. 

This implicit equation defines the curve shown in figure |5.3| . The figure does not show the 
isolated solutions 2, | ± ^z, — | ± ^i, which however are not in the spectrum, as they 



„8 „,2 



6 „A 



„4 „.6 



also come from solutions of (|5.3a| , |5.3b|) which satisfy ( |5.3c| ). 
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-1.5 -1 -0.5 0.5 1 
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Figure 2. Spectrum of u + u* + i(v + v* 



The intersection with the x-axis are the solutions at y — 0: 

x 16 - 6 x 15 + 5 x 14 + 32 x 13 - 59 x 12 - 40 x 11 + 136 x 10 - 32 x 9 - 48 x 8 

= (x 4 + 2x 3 -3x 2 -8x-3)(x-2) 4 x 8 = 

in particular, the spectral radius is the positive solution of 



x 4 + 2x 3 -3x 2 -8x-3 = 



which is 



1 



p(u 2 + u 3 ) = -{ : W3 (9 + (lSS-Gv^ 



/ 



(5.5) 



+ 3 



lSS-Gv 7 ^) 5 (135 + 6 v/249) ^ 



+ 



1.97148 



9 + (135 - 6 V249) 3 + (135 + 6 V249) 3 



With the exception of the double point at the apparent singularities of the curve only 
appear as such and the curve is actually smooth; for example, the curvature at the point of 
maximal modulus calculated in ( |5.5|) is approximately 166053.0, so the radius of curvature 
is about 1/166053.0 ~ 6.02219 x 10" 6 . 
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Example 5.2 (u + u* + i(v + v*)). Let u and v be freely independent Haar unitaries 
(e.g., u = X(gi), v = A (#2) , where g\ and gi are the generators of the free group F 2 on two 
generators. The distribution of the selfadjoint operator u + u* is the arcsine distribution 
with Cauchy transform G(() = — , 1 4 . The operator w + w* + i(v + t>*) recently came 



up as a candidate for a counterexample to the invariant subspace conjecture in the von 
Neumann setting; however, its "cousin", the circular element turned out to have plenty 
of invariant subspaces |pHO0 |. Let us compute the border of its spectrum here. The 
moment generating function of u + u* is 



s s 



VI - 4s 2 



from this we easily obtain the moment generating function of i(v + v* 

1 



git) = fiit) 



Now the L 2 -norm of the resolvent of a selfadjoint operator a is 

2 



1 — sa 



,-HI 2 



a 



,-1 



1 G(i)-G(§) 



and we have the equations 



sf(s)-sf(s) 

f(s) + sf'(s) 



Xsf-f-9 + 1 
sf-tg = 

/ 2 (1-4S 2 )-1: 



^ 2 (l + 4t 2 








Ims ^ 
Im s = 



A Grobner basis with respect to degrevlex order on (s, t, f, g) has 7 elements and the 
quotient algebra C[s,t, f,g]/I is four dimensional with basis ([g], [/], [t], [1]). Elimination 
did not succeed in this example and we chose the eigenvalue approach of section |5.1| . With 
the aid of the Grobner basis the matrix of multiplication by [s] is computed as 



M. 



Az 2 - 48 



z 5 + 24z 3 + 16z 
^ 4 - 24z 2 - 48 



z< 



+ 24z 3 + 16 



z 



2Az 2 + 16 



z A + 24^ 2 + 16 
z A + 16z 2 + 48 

z 5 + 24z 3 + 16z 



z 
1 

z 




1 

z 



z 2 + 16 



z 4 + 2Az 2 + 16 
2z 2 -8 

~ z 4 + 24z 2 + 16 
3z s + 12z 

z 4 + 2Az 2 + 16 
z 2 + 12 

~ z 4 + 2Az 2 + 16 



z 6 + 16z 4 - 48^ 2 + 128 
.+ 20z 5 - 80z 3 - QAz 



z 7 

z 6 - 16z 4 - A8z 2 - 12£ 

"2 7 + 20z 5 -802 3 -64z 
z 6 - 20^ 4 - 8O2 2 + 64 

z 6 + 20z 4 - 80z 2 - 64 
8z 2 + 32 

z 5 + 24^ 3 + 16z 



The result is shown in figure 5.3 
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Figure 3. Spectrum of U2 + u% + 

Example 5.3. As an example with more than two summands, let us consider U2+U3+U4. 
Generally, the equations for a sum u k, where Uk are free unitaries with = 1, are as 
follows. If, for given A, there is a solution Sj for the system 

1 



A 



1 _ n k 
1 b k 



I |2 1 I 4 . I |2n fc -2 

\sk\ + \sk\ H h \sk\ 



"fel 2 1 I I 2 1 I I 4 1 
- s k \ +\ s k\ +\Sk\ +■ 



+ \Sk 



2n h 



< 1 



then A is not in the spectrum. The example with n\ 
fig. 5.3. 



2, rti = 3, 113 = 4 is shown in 
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